Bose Condensate in a Double- Well Trap: 
Ground State and Elementary Excitations 
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We study the Bose-Einstein condensate in the MIT double- well trap. We calculate the ground 
state density profile of ^''Na atoms and the Bogoliubov spectrum of the elementary excitations 
as function of the strength of the double-well barrier. In particular, we analyze the behavior of 
quantum-mechanical collective excitations. Finally, we discuss the observability criteria for macro- 
scopic quantum tunneling and macroscopic quantum self-trapping. 
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Three years ago the Bose-Einstein condensation of al- 
kali vapors ^^Rb, ^'^Na and ''Li has been achieved in mag- 
netic harmonic traps at temperature of the order of 100 
nK [1-3]. Theoretical studies of the Bose-Einstein con- 
densate (BEC) in harmonic traps have been performed 
for the ground state [4-8], collective low-energy surface 
excitations [9-11] and vortex states [6,12]. 

In this paper we study the BEC in the MIT double- 
well trap given by a harmonic anisotropic potential plus a 
Gaussian barrier along the z axis, which models the effect 
of a laser beam perpendicular to the long axis of the con- 
densate. In the MIT experiment [13], the macroscopic 
interference of two Bose condensates released from the 
double minimum potential has been demonstrated. Such 
phenomenon has been theoretically reproduced [14] by 
using the Gross-Pitaevskii (GP) equation [15]. Here, we 
concentrate on the ground state properties of the con- 
densate and calculate the spectrum of the Bogoliubov el- 
ementary excitations as a function of the intensity of the 
laser field. A comparison between our calculations and 
future experiments will clarify the accuracy of the GP 
equation and the role of correlations in Bose condensates 
with up to 4 X 10^ particles. 

By varying the strength of the barrier one can ob- 
serve macroscopic quantum effects, like the formation of 
two Bose condensates, the collective oscillations and the 
quantum tunneling. 

The Gross-Pitaevskii energy functional [15] of the BEC 
reads: 

I = / rf'r l^jVM'Wl^ + Fe.tWlM'Wl^ + ^lvl/(r)l4 , 

(1) 

where ^'(r) is the wave function of the condensate nor- 
malized to unity, Vexti'^) is the external potential of the 
trap, and the interatomic potential is represented by a 
local pseudopotential so that g = 47r?i^ a^/m is the scat- 
tering amplitude (a^ is the s-wave scattering length) . N 



is the number of bosons of the condensate and m is the 
atomic mass. The extremum condition for the energy 
functional gives the GP equation 



- 2;^V2 + V,xt{r) + 5^1*(r)l2j ^,{v) = ^x^^{v) , (2) 

where /i is the chemical potential. This equation has the 
form of a nonlinear stationary Schrodinger equation. 

We study the BEC in an external potential with cylin- 
drical symmetry, which is given by 



^z2 + t/oexp(^) , (3) 



where p — \J + y^, z and the angle 9 are the cylindrical 
coordinates. The parameter values appropriate for Ref. 
[13] are = 27r x 250 Hz, = 27r x 19 Hz, and a = 
6 p,m. The anisotropic harmonic trap implies a cigar- 
shaped condensate (A = ujz/^jJp — 15/250 < 1), where z 
is the long axis, and the Gaussian barrier of strength [/q 
creates a double-well potential. 

We perform the numerical minimization of the GP 
functional by using the steepest descent method [18]. 
It consists of projecting onto the minimum of the func- 
tional an initial trial state by propagating it in imaginary 
time. At each time step the matrix elements entering the 
Hamiltonian are evaluated by means of finite-difference 
approximants using a grid of 200 x 800 points. In our 
calculations we use the z-harmonic oscillator units. For 
^^Na atoms, the harmonic length is az = {%/ {mLOz))^^"^ = 
4.63 //m and the energy is Tiujz — 0.78 peV. Moreover, we 
use the following value for the scattering length: Os = 3 
nm [3]. Most of our computations has been performed 
for iV = 5 X 10^ atoms, a value typical of the MIT exper- 
iment [13]. 

In Figure 1 we show the ground state density pro- 
file of the ^•^Na condensate for different values of the 
strength of the barrier. By increasing the strength, the 
fraction of ^^Na atoms decreases in the central region 
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and the Bosc condensate separates in two condensates. 
As shown in Table 1, the condensate shghtly expands in 
the z direction due to the barrier potential at the origin. 
The numerically calculated density profiles arc in good 
agreement with the phase-contrast images of the MIT 
experiment [13] and with the Thomas-Fermi (TF) ap- 
proximation, which neglects the kinetic term in the GP 
equation. Due to the large number of atoms involved 
(A'' = 5 X 10^), only near the borders of the wave func- 
tion there are small deviations from the TF approxima- 
tion. Note that the potential barrier Uq can be written 
as C/c/fcs = {^TnK)P/a'^ (/xm^/mW), where P is the 
total power of the laser beam perpendicular to the long 
axis of the condensate and a = 6/xm is the beam radius 
[19]. The conversion factor is 0.09 mW/(7ia;z), such that 
Uq = 100 (in fiLOz units) gives a laser Power P = 9 mW. 

Another important property of the BEG is the spec- 
trum of elementary excitations. To calculate the energy 
and wavefunction of the elementary excitations, one must 
solve the so-called Bogoliubov-de Gennes (BdG) equa- 
tions [20-22]. The BdG equations can be obtained from 
the linearized time-dependent GP equation. Namely, one 
can look for q angular momentum solutions of the form 

*(r,t) = e-^'^*[V'(/9,z)+e^'^'(u(p,z)e-*'"*+z;*(p,z)e'"*)] 

corresponding to small oscillations of the wavefunction 
around the ground state solution ip. By keeping terms 
linear in the complex functions u and v, one finds the 
following BdG equations 

Heffu{p, z) + gN\ip{p, z)fv{p, z) = Tiw u{p, z) , 

Heffv{p, z) + gN\il){p, z)fu{p, z) = -Tiw v{p, z) , (4) 

where 

ld_ 

+^^+Vea=t{p,z)-p + 2gN\i^{p,z)\'' . 

The BdG equations allow one to calculate the eigenfre- 
quencies tu and hence the energies huj of the elementary 
excitations. This procedure is equivalent to the diago- 
nalization of the N-body Hamiltonian of the system in 
the Bogoliubov approximation [22]. The excitations can 
be classified according to their parity with respect to the 
symmetry z —z. 

We have solved the two BdG eigenvalue equations by 
finite-difference discretization with a lattice of 40 x 40 
points in the (p, z) plane. In this way, the eigenvalue 
problem reduces to the diagonalization of a 3200 x 3200 
real matrix. We have tested our program in simple mod- 
els by comparing numerical results with analytical so- 
lutions and verified that a 40 x 40 mesh already gives 
reliable results for the lowest part of the spectrum [12]. 



In Table 2 we show the q ^ lowest elementary exci- 
tations of the Bogoliubov spectrum for the ground state 
of the system. When the Gaussian barrier is switched 
off, one observes the presence of an odd excitation at en- 
ergy quite close to hoj = 1 (in units hojz). This mode 
is related to the oscillation of the center of mass of the 
condensate along the z-axis, due to the harmonic confine- 
ment. This collective oscillation is an exact eigenmode of 
the problem characterized by the frequency w^, indepen- 
dently of the strength of the interaction. The inclusion of 
the Gaussian barrier modifies the harmonic confinement 
along the ^;-axis and this odd collective mode decreases 
by increasing the strength of the barrier (see Figure 2). 

In cylindrical coordinates, the other collective mode 
of the center of mass, due to the harmonic confinement 
along the radial axis, is an off- axial oscillation with an- 
gular quantum number q = 1. Such oscillation is the 
lowest even mode for = 1 and in the absence of a gaus- 
sian barrier, it is exactly equal to the radial frequency 
u)p {hujp = 250/19 = 13.158, in hujz units). One expects 
that this mode is only weakly affected by the gaussian 
barrier along the z-axis. In Table 3 are reported the first 
elementary excitations for q = 1. The lowest g = 1 exci- 
tation {Uw = 13.132) remains constant and differs by less 
than 2.5 °/oo from the theoretical prediction. Also when 
the BEG separates in two condensates, each condensate 
has the same off-axial (g' = 1) collective oscillation of the 
center of mass. 

As shown both in Table 2 and Table 3, for large values 
of the Gaussian barrier, i.e. when the BEG separates in 
two condensates, we find quasi-degenerate pairs of ele- 
mentary excitations (even-odd). The lowest q = mode 
and the ground-state of the GP equation constitute one 
of such pairs and get closer and closer as the barrier is 
increased. This is not surprising because in the infinite 
barrier limit we have two equal and independent Bose 
condensates with the same energy spectrum. 

An interesting aspect of BEC in double-well traps is 
the possibility to detect the macroscopic quantum tun- 
neling (MQT). The MQT has been recently investigated 
by Smerzi et al. [23,24]. They have found that the time- 
dependent behavior of the condensate in the tunneling 
energy range can be described by the two-mode equa- 
tions 

z = — z^ sm(j) , (j) = Az-\ — cos (/> , (5) 

Vl — -2 

where z = {Ni — N2)/N is the fractional population im- 
balance of the condensate in the two wells, </> = 0i — ^2 
is the relative phase (which can be initially zero), and 
A = 4:E^"*/AE°. E'""-* is the interaction energy of 
the condensate and AE'" is the kinetic -|-potential en- 
ergy splitting between the ground state and the quasi- 
degenerate odd first excited state of the GP equation. For 
a fixed A (A > 2), there exists a critical Zc = 2v'A — 1/A 
such that for < 2; << there are Josephson-like os- 
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dilations of the condensate with period t = tq/ a/I + A, 
where tq = 2iTh/A.E^. But for < z < 1 there is 
macroscopic quantum self-trapping (MQST) of the con- 
densate: even if the populations of the two wells arc ini- 
tially set in an asymmetric state {z ^ 0) they maintain 
the original population imbalance without transferring 
particles through the barrier as expected for a free Bosc 
gas. This two-mode approximation seems quite reliable. 
In fact, we have compared the predictions of the two- 
mode equations with the numerical solutions of the ID 
time-dependent GP equation in different regimes, finding 
a very good agreement (relative difference in the period 
of the Joscphson-like oscillations less than 1 °/o)- By us- 
ing the ID time-dependent GP equation, we have studied 
the dynamics of the condensate also outside the tunnel- 
ing region, i.e. when the chemical potential is higher 
than the gaussian barrier and the two-mode equations 
do not hold. In such case, starting, for example, with the 
condensate in one well, it docs not oscillate nor remains 
self-trapped but instead spreads over the two wells. 

By solving the stationary GP equation in the MIT 
double-well trap with ^"^Na, we find that the parameter 
A is larger than 10"* also when few particles are present. 
Nevertheless, we can control the dynamics of the conden- 
sate by reducing the scattering length ag and the thick- 
ness (7 of the laser beam. In Table 4 it is shown that, as 
expected, the parameter A scales linearly with a^. This 
is an important point because recently it was confirmed 
experimentally the fact that it is now possible to con- 
trol the two-body scattering length by placing atoms in 
an external field [25]. This fact opens the way to a di- 
rect observation of a macroscopic quantum tunneling of 
thousands of atoms through a potential barrier. 

Note that also the geometry and the dimensions of the 
trap play an important role. In fact, if the condensate is 
less cigar-shaped (greater A = LOz/iOp) the system has a 
lower chemical potential and weaker gaussian barriers are 
requested to enter the tunneling regime. Moreover, be- 
cause the strength of the nonlinear self-interaction scales 
as flg/a^, by increasing the dimensions of the trap we re- 
duce the effect of the nonlinearity thereby favoring quan- 
tum tunneling. 

In conclusion, we have shown that one can observe 
interesting macroscopic quantum-mechanical effects by 
studying the Bose-Einstein condensate in a double-well 
trap. We have accurately reproduced the formation of 
two Bose condensates observed by phase-contrast im- 
ages at the MIT experiment [13]. Moreover, by using 
the Bogoliubov-de Gennes equations, we have analyzed 
the behavior of elementary excitations as a function of 
the strength of the double-well barrier. In particular, 
we have identified two collective excitations which have 
a different fate by increasing the strength of the barrier: 
an odd q = {) mode that asymptotically goes to zero and 
a even g = 1 mode that remains constant also when the 
BEG separates in two condensates. We hope that our 



calculations will stimulate new, precise measurments of 
the excitation spectrum in this system in order to assess 
the validity of the theoretical framework adopted in the 
theoretical analysis. We have also considered the macro- 
scopic quantum tunneling. Our calculations suggest that, 
in the tunneling region and with ^^Na atoms, one sees 
only the macroscopic quantum self-trapping (MQST) of 
the condensate also when a small laser-sheet thickness is 
applied. We have shown that to get outside the MQST 
regime it is necessary to strongly reduce the scattering 
length or to increase the dimensions of the trap. Note 
that at nonzero temperature, BEG depletion and ther- 
mal fluctuations will slightly modify the parameters of 
the tunneling and will damp the coherent oscillations. 
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FIG. 2. Lowest elementary excitation hu)[ ' with q = vs 




FIG. 1. Particle probability density in the ground state 

of = 5 X 10^ ■^^Na atoms as a function of the z axis at 
r = (symmetry plane) . The curves correspond to increasing 
values of the strength Uo of the barrier (from to 500), in 
units of hujz = 0.78 peV. The laser power is given by the 
conversion formula P = 0.09 x Uo mW. Lengths are in units 
of Oz = 4.63 /Ltm. 



barrier energy Uo for 
Fig. 1. 



5 X 10 Na atoms. Units as in 



Up E/N ^ V<p'^> V< > 



211.142 294.909 
100 221.334 303.910 
200 227.467 310.705 
300 230.852 315.297 
400 232.855 317.952 
500 234.283 319.763 



0.989 


9.145 


0.982 


9.635 


0.986 


10.011 


0.993 


10.268 


0.997 


10.417 


0.999 


10.518 



Table 1. Ground state of iV = 5 x 10^ '^^Na atoms. 
Energies are in units of hojz = 0.78 peV (w^ = 27r x 19 
Hz). Lengths are in units of Uz = 4.63 ^m. The laser 
power is given by the conversion formula P = 0.09 x Uq 
mW. 
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1.000 


1.580 


2.120 


2.645 


3.164 


100 


0.939 


1.596 


2.022 


2.640 


3.004 


200 


0.799 


1.624 


1.886 


2.711 


2.916 


300 


0.311 


1.643 


1.672 


2.724 


2.744 


400 


0.003 


1.655 


1.655 


2.730 


2.730 


500 


10-^ 


1.663 


1.663 


2.744 


2.744 



Table 2. Lowest elementary excitations of the q = 
BdG spectrum for the ground state of A'' = 5 x 10^ 
^■'Na atoms. and mean odd and even z-parity, 
respectively. Units as in Tab. 1. 



i+r 



FT 



13.132 
100 13.132 
200 13.132 
300 13.132 
400 13.132 
500 13.132 



13.165 
13.158 
13.145 
13.133 
13.132 
13.132 



13.214 
13.218 
13.222 
13.225 
13.226 
13.227 



13.278 
13.260 
13.236 
13.225 
13.226 
13.227 




200 300 
barrier energy 



500 



Table 3. Lowest elementary excitations of the q = 
1 BdG spectrum for the ground state of A'' = 5 x 10^ 
^^Na atoms. and mean even and odd z-parity, 
respectively. Units as in Tab. 1. 



as/ a 



To (sec) Zc 



10"^ 1108.337 14.583 0.060 

10-2 133.643 13.887 0.173 

10-3 1.390 13.842 none 

10-^ 0.103 10.253 none 



Table 4. Parameters of the MQT for different values 

of the scattering length with a^" = 3 nm and tq = 
27r?l/Ai^°. Condensate with N = bx 10^ atoms. Barrier 
with Uq = 20 and cr = 1.5 /xm. Unit of Uo as in Tab. 1. 
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